Journal of Molecular Graphics and Modelling

The NAD+ dependent deacetylase sirtuin-1 (SIRT1) is known to elicit cellular defenses against aging, cancer, and other aberrant pathologies. Previous studies have identified an intrinsically disordered region of SIRT1 comprised of N-terminal residues 1-52, herein referred to as motif A, which activates SIRT1 activity, likely through intramolecular interactions. Additionally, phosphorylation of N-terminal residues Ser27 and Ser47 has been shown to be important for regulating SIRT1 activity and stability. The lack of in vitro characterization of these effects hampers our further understanding of the role of motif A in SIRT1 regulation. In this study, we elucidate the role phosphorylation plays in motif As structure as well as its regulatory effects on SIRT1 activity against Ac-p65. We find that phosphomimetic mutation at Ser27 significantly increases the activation effect of motif A towards SIRT1. This result is supported by molecular dynamics simulations of the phosphomimetics, which reveal stabilization of different transient structures for motif A depending on whether Ser27 and Ser47 have been modified. A key finding suggested by this study is that phosphorylation of S27 appears to activate SIRT1 by causing motif A, which is intrinsically disordered in the WT, to fold into an ordered structure. This conclusion is based on both the experimental findings and simulation results. These findings contribute to our understanding of SIRT1 regulation, specifically the role played by phosphorylation within the N-terminal disordered region.

The Na+/K+-ATPase (NKA) regulates ion balance in the kidney and influences cellular processes like proliferation and apoptosis through its signal transduction. The endogenous ligand 20-Hydroxyeicosatetraenoic acid (20-HETE) contributes to inflammation and fibrosis in chronic kidney disease (CKD) and inhibits NKA activity in renal tubules. However, the molecular mechanism of this interaction remains unclear. In this study, we used in-silico approach to investigate the potential interaction between 20-HETE and NKA. Various ligands, including known NKA ligands such as cardiotonic steroids (CTS), 20-HETE, and negative controls, were docked using rigid and Induced Fit Docking to predict the affinity of the ligands toward NKA. Binding free energy calculations with the Prime Molecular mechanics with generalized Born and surface area (Prime MM/GBSA) tools were used to confirm the involvement of key amino acids in ligand-receptor interactions. The docking analyses revealed that 20-HETE exhibited a binding affinity comparable to negative control, with some differences between rigid and induced fit docking. Binding free energy data highlighted key amino acids in the 20-HETE and NKA interaction. Interaction fingerprint and mutations such as Ala330Gly and Val329Ala significantly reduced binding free energy, while Thr804Ala showed a notable decrease, underscoring the potential importance of these amino acids in ligand stabilization. These findings provide computational evidence supporting potential direct interaction between 20-HETE and NKA and identify candidate residues for future experimental validation.

The M. tuberculosis ESX-1 secretion system EccA1 enzyme is involved in the secretion of virulence factors and is essential for virulence and bacterial survival within the phagosome. Development of the small molecular inhibitors abolishing EccA1 function can yield new antivirulence drugs. In this study, we modeled the full-length EccA1 (573 residues, Mw [~]62.4 kDa) structure, which contains N-terminal TPR domain and a C-terminal CbxX/CfqX type ATPase domain. We have identified five ZINC compounds having binding energy i. e. Z1 (ZINC000004513760, -43.45 kcal/mol), Z2 (ZINC000000001793, -49.56 kcal/mol), Z3 (ZINC000005390388, -55.83 kcal/mol), Z4 (ZINC000257294577, -52.33 kcal/mol), Z5 (ZINC000004824264, -44.44 kcal/mol) through virtual screening of the ZINC compounds targeting C-terminal ATPase pocket of EccA1. The Z1-Z5 compounds were compared with ADP substrate having binding energy (Adenosine diphosphate, -35.00 kcal/mol), p97 ATPase inhibitors i.e. NMS873 (3-[3-cyclopentylsulfanyl-5-[[3-methyl-4-(4 methylsulfonylphenyl)phenoxy]methyl]-1,2,4-triazol-4-yl]pyridine, -48.68 kcal/mol), and CB5083 (1-[4-(benzylamino)-5H,7H,8H-pyrano[4,3-d]pyrimidin-2-yl]-2-methyl-1H-indole-4-carboxamide, -50.88 kcal/mol) against EccA1. The Z1-Z5 compounds exhibited good Absorption, Distribution, Metabolism, and/or Excretion properties (ADMTE). Pharmacokinetic properties and Lipinskys rule of five for Z1-Z5 compounds showed drug-like properties. 100 ns dynamics simulation analysis on EccA1 complexed with (i) Z1-Z5 compounds (ii) ADP substrate and (iii) NMS873 and CB5083 inhibitors showed high stability and biologically relevant conformation during dynamics simulation. These data indicate that Z1-Z5 compounds may act as potential inhibitors against EccA1 and provide avenues for new antivirulence drug development after in vitro and in vivo clinical trials.

Formation of the {beta}-amyloid (A{beta}) plaques is a pathological hallmark of Alzheimers disease (AD), and is believed to be a primary cause of dementia in elderly individuals. In the present work, we simulated the conformational evolution of A{beta}42 dimers in solution and in membrane-like environment to explore the folding of A{beta}42 along fibrillation. The molecular dynamics (MD) simulation was steered by experimental internuclear distance restraints obtained using solid-state nuclear magnetic resonance (ssNMR) spectroscopy. Our results revealed that several hydrophobic and polar motifs within the A{beta}42 sequence played key roles in the early-stage nucleation process of fibrillation and those motifs are also the stabilizing agents in the mature fibrils judged by the energy contribution. Our results also indicated that the peptide association with membrane bilayers could modulate the structural evolution pathways towards fibrillation. These findings contributed to a better understanding of the molecular level structural polymorphisms inherent to A{beta}42 fibrils. Further, the current work demonstrated that the combination of MD simulations with ssNMR-based experimental restraints provided a reliable method for studying structural changes of A{beta}. HighlightO_LIUsing solid-state NMR restraints guided molecular dynamic simulation, {beta}-amyloid dimers displayed consistent {beta}-strand-prone regions, which are major stabilizing segments for mature fibrils. C_LIO_LI{beta}-amyloid dimers evolved differently with or without interacting with the lipid bilayers. C_LIO_LIExperimental restraints guided simulation provided molecular level insights about early-stage interactions along the progress of {beta}-amyloid fibrillation C_LI

Rademikibart (CBP-201) is a human monoclonal antibody with higher binding affinity to IL-4R compared to dupilumab. Dupilumab is a first-generation interleukin-4 receptor alpha (IL-4R) inhibitor for treating IL-4R-dependent inflammatory disorders, including several dermatologic and respiratory conditions. Rademikibart, however, demonstrated better inhibition of STAT6 intracellular signaling in vitro and similar potency in inhibiting both IL-4 induced TARC release and IL-4 induced B cell activation. To further characterize the molecular function of rademikibart and its differentiation from dupilumab, we determined the crystal structure of the rademikibart fragment antigen binding (Fab) bound to IL-4R at 2.71 [A] resolution and compared this to the 2.82 [A] resolution structure of dupilumab Fab bound to IL-4R. The rotation angle between dupilumab and rademikibart bound to IL-4R is 54.88{degrees}. This rotation enables the binding epitopes of rademikibart, but not dupilumab, on IL-4R to overlap more closely with the conserved binding interface naturally utilized by IL-4 and IL-13 cytokines. Molecular dynamics (MD) studies on rademikibart and dupilumab bound to IL-4R examined the stability of the complexes and effects of amino acid mutations on receptor complex formation. MD simulations demonstrated that the third interface loop (residues 145 to 153 in domain 2) of IL-4R interacts directly with rademikibart, which is absent in the dupilumab/IL-4R complex. This finding is confirmed by increased hydrogen bond interactions at the interface between rademikibart and IL-4R, demonstrating superior binding energy for rademikibart. Through analysis of the x-ray crystallography structures, MD-equilibrated structures, and computational point-mutation analysis of rademikibart, we identified residue Y50 and R55 of the light chain and R97, R99, and Y101 of the heavy chain of rademikibart as key residues interacting with IL-4Rs third interface loop. Our data provides a molecular and structural rationale for the enhanced IL-4R inhibition by rademikibart over dupilumab, confirming rademikibart as an optimized second-generation IL-4R inhibitor.

Marburg virus (MARV) is a highly pathogenic filovirus that causes hemorrhagic fever with a high mortality rate, with very limited treatment options. The urgent need for targeted antiviral agents emphasizes the importance of structure-based drug discovery approaches. The present study aimed to evaluate the antiviral potential of Withaferin A (PubChem CID-265237) against three key proteins of MARV: viral protein 35 (VP35), and nucleoproteins (NP). Three-dimensional structures of these proteins were retrieved from RCSB-Protein Data Bank and docked with Withaferin A using AutoDock Vina. The ligand demonstrated favourable binding affinities towards all three viral targets, indicating strong interaction potential at functionally relevant sites. Drug-likeness and pharmacokinetic properties predicted using SwissADME and pkCSM indicated acceptable ADMET profiles that comply with key drug-like criteria. To validate the stability of the docking, molecular dynamics simulations (GROMACS, 100 nanoseconds) were conducted. The protein-ligand complexes exhibited stable root mean square deviation (RMSD), root mean square fluctuation (RMSF), and consistent hydrogen bonding patterns throughout the simulation. The MM-GBSA binding free energy analysis further supported favorable binding energetics, predominantly driven by van der Waals and electrostatic interactions. Altogether, these findings demonstrate that Withaferin A exhibits promising multi-target inhibitory potential against key MARV proteins. This study provides molecular insights into ligand-protein interactions and supports further experimental validation of Withaferin A as a potential therapeutic candidate against Marburg virus.

Signaling through the insulin receptor (IR) and the type 1 insulin-like growth factor receptor (IGF1R) is modulated by secreted hormones and growth factor ligands (e.g. insulin and insulin-like growth factor 1, IGF1). Impaired signaling in these receptors often leads to diabetes and oncogenic diseases. The discovery of entirely novel viral insulin/IGF-like peptides (VILPs) that can stimulate receptors from the insulin family has raised questions about their structures and binding modes to receptors. These peptides exist in a single-chain (sc) or a double-chain (dc) configuration with folds likely similar to IGF1 and insulin, respectively. The interactions of VILPs with the human receptors are beginning to be mapped but little is known about their interactions with the receptors in fish-the host organism for viruses known to carry these peptide sequences. We have previously reported [Chuard et al., Cell Rep. 2025 44(8):116149] structural models of several VILPs from the Iridoviridae virus family bound to their cognate receptors in Zebrafish (Zeb). In this work, we conducted all-atom molecular dynamics (MD) simulations of these peptides and their receptor-bound complexes along with free energy calculations to assess the energetic contributions of VILP residues for their binding to Zebrafish receptors. Most of the observed Zeb insulin/Zeb {micro}IR and Zeb IGF1/Zeb {micro}IGF1R site 1 interactions are consistent with previously known interactions of human peptides with their receptors, highlighting similarities in their binding modes. However, we also report some non-conserved residues in VILPs that establish significant and unique interactions with residues in Zeb receptors. Furthermore, we identified residues in each VILP which can be potentially mutated into conserved insulin/IGF1 residues to possibly enhance the binding affinity of these peptides.

Like other proteins, monoclonal antibodies - important biodrugs- are subject to post translational modifications, especially the N-glycosylations. However, the effect of the N-glycosylations remains poorly studied and atomistic details about their influence are rarely available. Moreover, the few existing studies focus on the prevalent immunoglobulin G1. To go further in the understanding of the impact of glycosylations, we have carried out a comparative exploration of the effect of N-glycosylations on two different classes of antibodies, namely Mab231, an IgG2 and the pembrolizumab, an IgG4. The two antibodies differ by their sequences, their length, their 3D structure but also by the location and composition of the glycans. In the present work, detailed and important information were gained through molecular dynamics simulations where both monoclonal antibodies were studied without and with the presence of their glycans. The results of 1.5 {micro}s of sampling for each system show that glycosylation does not drastically alter the overall conformational landscape of either antibody, whatever the metrics considered. However, it measurably modulates local flexibility, inter-domain correlated motions, and the relative orientation of the Fab arms with respect to the Fc domain, with statistically significant shifts in key geometric descriptors. Importantly, contact analysis reveals that glycan interactions extend beyond the Fc region to reach Fab residues. The allosteric network calculations demonstrate that the influence of Fc-bound glycans propagates even until the Fab framework regions in both mAbs, which could impact the antigen binding. The nature and magnitude of these effects are subclass-dependent, reflecting differences in glycan composition, hinge architecture, and three-dimensional organization Our findings challenge the prevailing view that Fc glycosylation uniformly promotes CH2 domain opening. More importantly, it underscores the necessity of considering full-length structures and IgG subclass diversity in glyco-engineering strategies.

Toll-like receptors (TLRs) are vital components of the innate immune system, recognizing both exogenous pathogens signals (PAMPs) and internal stress signals (DAMPs). TLR2 is unique among the human (Homo sapiens) TLR family members, as it contains a large cavity for binding hydrophobic ligands, such as lipoteichoic acid (LTA) and di/triacyl lipopeptides (Pam2/3CSK4). This study analyzed the structural phylogeny of cavity presence in the TLR2 lineage in vertebrates (vTLR) enabled by AI protein structure predictions and explored the potential convergent evolution of similar features in invertebrates (iTLRs). Analysis of AI models of TLR2s shows that this cavity is consistently present in TRL2 orthologs across jawed vertebrates (Gnathostomata). In jawless vertebrates (Cyclostomatha), these cavities were found in lamprey (Petromyzon marinus) TLR2 model, but only in some extant hagfish (Myxini), suggesting an ancestral origin in basal vertebrates followed by lineage-specific losses. TLR2 paralogs were found in several species, with a similar central cavity but potentially different ligand specificities. In silico ligand docking showed Pam2CSK4 binds to this cavity in all TLRs and paralogs consistently, demonstrating the conserved function of the ligand-binding pocket in gram-positive bacteria recognition across TLR2 branches. Changes in the TLR2 cavity size and shape in some vertebrate groups show the evolution of this DAMP recognition mechanism adapted to its respective pathogens. iTLRs form a separate phylogenetic branch with distinct structural features, but in literature some are considered to be TLR2 orthologs. Indeed, TLRs from some species of Helobdella and Ciona, contain a cavity with some similarity to that in the vTLR2 lineage. However, detailed structural comparisons of their location in the LRR domain and the structural details of the models suggest that their cavities have developed independently from that in TLR2s. Smaller cavities are present in other branches of the LRR family, but show different locations, shapes, and features, indicating that the binding of small ligands in the internal cavities within the LRR domains evolved multiple times in the LRR domain family history.

Pseudomonas aeruginosa, a Gram-negative opportunistic pathogen is well known for life-threatening acute infections among the human population. The bacterium can withstand most antibiotics by using their high levels of inherent and acquired resistance mechanisms such as Biofilm-EPS, Persistence, and Quorum sensing (QS). Owing to the importance of adaptive antibiotic multi-drug resistance of P. aeruginosa, the current investigation is aimed to explore the phytochemicals derived from mangrove plants as potential agents to control biofilm and drug resistance mechanisms through a multi-mechanistic computational approach. For identifying potential compounds and target, In-silico drug repurposing technique is implemented by docking/virtual screening of 49 phytochemical compounds against 18 proteins involved in the Persister Cell formation, QS, and EPS synthesis in P. aeruginosa which resulted the proteins RelA and SpoT (persistence), PqsA, and PqSR (QS), and PelA and PelB (EPS synthesis) and compounds Taraxerone and Taraxerol to be potential. The results of docking were well corroborated with MD simulations. These targets and compounds explored through in-silico approach, are found to target potential antimicrobial pathways involving EPS synthesis, persistence genes, and QS, aiming to enhance antibiotic efficacy. Further, this study could be reference for in-vivo and in-vitro investigations to evaluate the further effectiveness of the compounds and potentiality of the proteins for MDR therapeutics of P. aeruginosa.

Flavonoids have been widely investigated for their antiviral and anti-inflammatory properties, but their mechanisms of action often remain insufficiently defined. In the present study, high-purity flavonoids were evaluated using an integrated workflow combining molecular docking, LigPlot+ interaction mapping, surface plasmon resonance (SPR), fluorescence-based TMPRSS2 inhibition assays, and cell-based viability studies. Docking with AutoDock Vina identified Hesperidin as the strongest overall candidate among the compounds evaluated. Hesperidin showed strong active-site engagement with TMPRSS2, including interactions with catalytic residues His296, Asp345, and Ser441, and stable binding within the SARS-CoV-2 main protease (Mpro) pocket. Comparative docking showed weaker or more peripheral interaction patterns for Rutin and moderate Spike binding for Hesperidin and Rutin. Experimental validation demonstrated dose-dependent inhibition of TMPRSS2 activity with an IC50 of 79.1 {micro}M for Hesperidin and 43.5 {micro}M for Hesperetin, while Rutin showed partial inhibition without a defined IC50 in the tested range. In Calu-3 cells, pre-treatment with Hesperidin or Rutin reduced SARS-CoV-2 Spike-induced cytotoxicity by approximately 30% without detectable intrinsic toxicity at the concentrations tested Docking analysis of Hesperidin and Rutin with the SARS-CoV-2 Spike protein revealed moderate interaction patterns involving residues such as Asn343, Ser371, and Val367. Hydrogen bond distances were generally in the range of approximately 2.9-3.3 [A], indicating moderate stabilization compared with the stronger active-site interactions observed for Hesperidin in TMPRSS2. The resulting binding poses suggest that these flavonoids can associate with structurally relevant regions of the Spike receptor-binding domain; however, they do not strongly overlap with the key residues required for ACE2 interaction. Rutin, in particular, exhibited a more peripheral and distributed binding mode within the Spike-ACE2 complex, indicating limited potential for direct disruption of the binding interface. In addition to SARS-CoV-2 targets, docking analysis extended to influenza viral proteins revealed moderate interaction of Hesperidin with hemagglutinin (HA) and strong catalytic-pocket binding of Rutin to neuraminidase (NA), involving key residues associated with enzymatic activity. These findings broaden the scope of the study to include influenza viral entry and release mechanisms, supporting a multi-virus, multi-target framework.

Synaptic vesicle glycoproteins 2 (SV2) are integral membrane proteins essential for neurotransmitter release and are implicated in neurological disorders including epilepsy and Parkinsons disease. In the attempt to develop a ligand selective for SV2C, and in collaboration with UCB, UCB-F was identified as a potential candidate. However, the affinity of UCB-F to SV2C was found to be temperature dependent, decreasing by about 10-fold from +4 to 37 degrees. UCB1A was subsequently identified as SV2C ligand displaying in vitro a 100-fold selectivity for SV2C compared with SV2A. In this study we investigated whether the binding of UCB-1A to SV2A and SV2C was affected by the temperature. A combination of experimental binding assay data and molecular dynamics (MD) simulations were used. The binding studies revealed that UCB1A affinity for SV2A decreased significantly at 37 {degrees}C compared with 4 {degrees}C, whereas binding to SV2C remained largely unchanged. MD simulations reproduced these observations, namely that ligand RMSD values at 310 K showed that UCB1A binding fluctuated markedly in the SV2A complex, with many trajectories exceeding the 3.0 [A] stability cutoff, whereas UCB1A remained relatively well-anchored in SV2C under the same conditions. Structural analysis showed that, while UCB1A adopts a conserved binding pose across all isoforms stabilized by {pi}- {pi} stacking and a hydrogen bond with Asp, SV2C possesses a unique stabilizing feature. In SV2C, Tyr298 is less exposed to the solvent and engages in a persistent hydrogen bond with Asparagine, a structural feature that reinforces pocket stability and limits temperature-induced destabilization. This interaction is absent in SV2A, consistent with its greater temperature sensitivity. Together, these findings provide a mechanistic explanation for the experimentally observed temperature independence of UCB1A binding to SV2C. More broadly, the results highlight the importance of incorporating physiologically relevant temperatures into SV2 ligand evaluation and demonstrate how combining experiments with simulations can uncover isoform-specific mechanisms of ligand recognition and stability.

Background/ObjectivesThe unprecedented structural and binding data for antibodies to the SARS-COV2 virus taken together with the mutations for the spike protein allows for a broad simulation study of antibody-spike protein binding. This provides an understanding of the co-evolution of human immunity and viral immunity escape. MethodsWe utilized the YASARA molecular dynamics program to generate initial antibody-spike structures and simulate to equilibration for six SARS-COV2 variants and 10 different antibodies sampling two different binding regions to the receptor binding domain of the spike (especially for the Class I antibodies in the same part of the spike which attaches to the ACE2 receptor protein) and one to the N-terminal of the spike. Starting structures for antibody binding to variant spike proteins are perturbatively achieved through point mutations and insertions/deletions in the YASARA program. We employed YASARA to measure interfacial hydrogen bound counts between antibodies and variant spike proteins, and the HawkDock MMGBSA program to characterize trends in binding energies with mutation for four of the antibodies. We utilized the VMD program to analyze the time course of hydrogen bond populations. ResultsAs seen in previous studies, interfacial hydrogen bond counts serve as an excellent proxy for binding energies without the large systematic error inherent in the latter. We find that there is generally a decline in antibody binding strength, as measured by interfacial hydrogen bond counts, with viral evolution, but that a modest re-entrance of binding strength is present for most antibodies studied. Generically, the antibody heavy chain binds more strongly to the spike protein, through for approximately half the antibodies the light chain binding strength converges to the heavy chain strength with viral evolution. ConclusionsThe key conclusion is that the identified re-entrant immunity, speculatively arising from a balancing of maintenance of ACE2-spike binding while escaping antibodies through mutation, allows for some maintenance and even strengthening of immunity for later viral strains from early infection or vaccination.

The intracellular transport system is pivotal for cellular function and integrity, facilitated by cytoskeletal motor proteins such as dynein, which traverse along microtubules (MTs). The heterogeneity of the tubulin isotypes composing MTs introduces functional diversity, potentially affecting cytoskeletal motor proteins interactions with the MT. This in silico study investigated the influence of amino acid sequence variations in the C-terminal tails (CTTs) of six different Homo sapiens tubulin isotypes, TUBB2A, TUBB2B, TUBB2C, TUBB3, TUBB4A, and TUBB5, highly expressed in human brain tumors, and assessed the isotypes effect on the binding of motor protein dynein to MT. Among these isotypes, TUBB2A, TUBB2B, and TUBB2C were found to affect conformational motions of the dyneins microtubule-binding domain (MTBD) and stalk domain. The investigation highlighted the novel role of isotype-specific variations in lateral interactions between tubulin protofilaments (PFs) in determining the proximity of the {beta}-CTT of the adjacent PF to the MTBD, potentially affecting dyneins motility and suggesting how changes in isotype expression directly influence dyneins velocity and processivity and contribute to transport defects associated with neurological disorders and cancers. Isolating specific tubulin isotypes experimentally is challenging due to their high sequence similarity and complex interactions with other microtubule-associated proteins. This makes it challenging to distinguish between different tubulin isotypes and their effects, particularly in tissues where multiple isotypes are co-expressed. Additionally, these isotypes are heavily modified in vivo by post-translational modifications, which further complicate the isolation of a single, unmodified tubulin isotype. As a result, computational approaches have been necessary in this study for exploring these effects in a controlled, isotype-specific manner.

We have conducted all atom molecular dynamics simulations of POPC and DPPC lipid bilayers using AMBER Lipid21 force field with eight different water models, including SPC/E, TIP3P, TIP3P-FB, TIP4P-FB, TIP4P-Ew, TIP4P/2005, TIP4P-D, and OPC, to identify the most compatible one without any modification. A number of parameters have been computed in order to understand the structure of the lipid bilayer: Area per lipid, Isothermal compressibility modulus, average Volume per lipid, electron density profile, bilayer thickness, X-ray and neutron scattering form factors, deuterium order parameter, and radial distribution function. The estimated Area per lipid, Isothermal compressibility factor, volume per lipid and bilayer thickness are highly consistent with experimental results for the SPC/E water model, indicating its suitability with the AMBER Lipid21 force field, insted of any modification. The bilayer electron density profiles of both the lipid bilayers demonstrate a little augmentation of water penetration with respect to the membrane surface for TIP4P-D water model. However, the experimental X-ray and neutron scattering form factors are aligning well with the simulated results for all studied water models, and TIP4P-D shows better for X-ray data. The deuterium order parameter for lipid acyl chains value less than 0.25 for all observed water models, depicting their disorderness for both the lipid bilayers. The lateral diffusion and reorientation autocorrelation function of the lipid molecules in both the bilayers are computed to reveal their dynamics across all water models. In comparison to other water models, the simulated trajectories predict better structure and reasonably fair dynamic properties for the SPC/E water model. The TIP4P-Ew water model reproduces the lateral diffusion co-efficient in close agreement with experiment. Reorientational dynamics for both the lipids in the bilayers for eight different water models are observed; the presence of slow and slowest time components corresponds to the lipid axial motion (wobble motion) and Twist/Splay motions. So, in view of the overall performance of the different water models with the AMBER Lipid21 all atom force field in reproducing membrane physical properties, the SPC/E water model appears to be an optimal choice.

The nuclear receptor superfamily is comprised of ligand-regulated transcription factors that contain an intrinsically disordered domain at the amino-terminal end, known as the N-terminal domain (NTD). While this poorly conserved domain is known to possess ligand-independent activation function (AF-1), few NTD functions are conserved between nuclear receptors (NRs). Identified roles in other receptors include androgen receptor (AR), estrogen receptor (ER) and mineralocorticoid receptor (MR). Here, we aim to define the function of the NTD of the farnesoid X receptor (FXR), a crucial regulator of lipid and bile acid metabolism. We show that the NTD engages in interdomain contact with other FXR domains. We also observe that the NTD interacts directly with coregulator proteins. Using mutagenesis, mammalian two-hybrid assays and molecular dynamics simulations, we identify and validate a novel SXXLF motif in the NTD which mediates interactions with both coregulators and the ligand binding domain. Mutation of the motif induces large changes in conformational and allosteric coupling in FXR. Our study identifies a new nuclear receptor-interacting motif that modulates the transcriptional activity of FXR. Graphical AbstractFXR-NTD regulates transcriptional activity through interdomain communication with the LBD and is also involved in co-activator recruitment. The SENLF motif is the first defined functional element within the FXR-NTD and mediates both NTD-LBD interaction and selective co-activator engagements to drive NTD-mediated transcriptional activity. O_FIG O_LINKSMALLFIG WIDTH=135 HEIGHT=200 SRC="FIGDIR/small/724725v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@5a37aorg.highwire.dtl.DTLVardef@2fa9e1org.highwire.dtl.DTLVardef@13a19daorg.highwire.dtl.DTLVardef@1775ed2_HPS_FORMAT_FIGEXP M_FIG C_FIG

BackgroundRising antimicrobial resistance rates, require new therapeutic approaches such as antimicrobial peptides (AMPs), which are part of the innate immune defense, as alternatives to antibiotics. In this study, we aim to unravel the antibacterial activity of human histone H1.2 peptide against Pseudomonas aeruginosa and its potential immune modulatory role. MethodsWe used a hemofiltrate peptide database for antimicrobial peptide prediction to identify novel human AMPs. Thirteen sequences of histone H1 were identified as putative AMPs, synthesized, and tested against bacterial ESKAPE pathogens in a radial diffusion assay. SYTOX green assay, electrophoretic mobility shift assay, and differential proteomics assays were conducted to determine the mode of action of H1.2 peptide fragment. A crystal violet assay was performed to evaluate the inhibition of biofilm formation. The cytotoxicity of the peptide was tested in LDH and Alamar assays. Finally, to visualize the contributions of H1.2 in NETs formation, scanning electron microscopy was performed. ResultsThe H1.2 peptide inhibited the growth of P. aeruginosa in a dose and pH-dependent manner without cytotoxicity towards mammalian THP-1 cells. It acts on intracellular targets to inhibit the growth of P. aeruginosa. STRING analysis from the differential proteomics assay showed that H1.2 targets the downregulation of proteins involved in the biogenesis of outer membrane proteins, including the folding and trafficking of outer membrane proteins across the cytoplasmic membrane. Scanning electron microscopy images showed that H1.2 forms NET-like structures capable of trapping and immobilizing P. aeruginosa. ConclusionThe characterized antimicrobial activity of H1.2 points to a role for human histone H1 fragments in innate immunity and may represent a promising approach for the development of novel antibacterial therapies. Graphical Summary O_FIG O_LINKSMALLFIG WIDTH=192 HEIGHT=200 SRC="FIGDIR/small/724237v1_ufig1.gif" ALT="Figure 1"> View larger version (36K): org.highwire.dtl.DTLVardef@1778ddborg.highwire.dtl.DTLVardef@26430org.highwire.dtl.DTLVardef@ffbfa2org.highwire.dtl.DTLVardef@7e38ae_HPS_FORMAT_FIGEXP M_FIG C_FIG Sec transport and BAM complex system including chaperone proteins and quality control proteases are inhibited by H1.2 in Pseudomonas aeruginosa.Outer membrane proteins (OMPs) are synthesized in the cytoplasm and transported across the inner membrane via the Sec translocase, assisted by SecA/SecB or ribosomes. In the periplasm, they are escorted by chaperones such as SurA to the BAM complex for insertion into the outer membrane. Here, we show that H1.2, an antimicrobial peptide, targets membrane biogenesis in P. aeruginosa through downregulating Sec translocase (SecA/SecB and SecYEG), SurA, and BAM complex. Therefore, leading to improper transfer, folding and insertion of OMPs into the outer membrane. Normally, misfolded proteins are degraded by the protease MucD to prevent toxic aggregation in the bacteria. However, with H1.2 inhibiting MucD the proteotoxic stress is exacerbated, ultimately compromising bacterial homeostasis and viability. Figure created using BioRender.com.

WD40 domains share a widespread {beta}-propeller fold, and often act as versatile scaffold proteins. Despite their central role in organizing dynamic cellular complexes, the molecular and structural mechanisms of many WD40 proteins remain poorly understood. Among them, DCAF7, an ubiquitously expressed and essential gene in human, also encodes a highly conserved WD40 protein in eukaryotic organisms. It is known to interact with multiple and functionnally diverse partners to coordinates cellular activity of several protein kinases as well as transcriptional regulators, thereby modulating key cellular processes such as cell growth, differentiation, and transcriptional regulation. However, the precise mode of action of DCAF7 is unknown and its important divergence in sequence from better characterize WD40 prevent information transfer by similarity. Structural interactomic can reveal how protein-protein interactions (PPIs) occur within an organism and are essential for understanding biological functions and developing new therapeutic strategies. Using SLiMAn2, AlphaFold2/3 and PSSMsearch, we identified a conserved -helical short linear motif (SLiM) in several well known DCAF7 partners that binds to the top surface of its {beta}-propeller. This motif was subsequently used to generate a regular expression, to identify potential new direct binders across the DCAF7 meta-interactome and the human proteome. Domain-domain interactions were also predicted for some other partners. Finally, modeling of oligomeric complexes with such new hits reveals the structural basis of DCAF7 scaffolding, with links to neurodevelopmental disorders such as autism.

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are are a family of voltage-gated, cyclic-nucleotide modulated Na+/K+ channels that regulate spontaneous rhythmic electrical activity in both the heart and the brain. Understanding differences in the responsiveness to cyclic adenosine monophosphate (cAMP) modulation between HCN isoforms would offer insight into the specific binding interactions that drive channel activation. Using all-atom molecular dynamics (MD) simulations and the free-energy perturbation (FEP) approach, we determined the absolute binding free energy of cAMP to the the cyclicnucleotide-binding domain (CNBD) of HCN isoforms 1-4. By studying the free-energy of ligand binding to the various isoforms of HCN, our study advances the understanding of HCN channel activation and modulation mechanisms. Overall, our work offers insight into explaining differences in channel sensitivity across the isoforms of HCN.

Chlamydia trachomatis is an obligate intracellular Gram-negative pathogen responsible for sexually transmitted infections and trachoma in humans. Although antibiotics are generally effective against acute infections, persistent chlamydial forms often exhibit reduced susceptibility during chronic infection. Chlamydia relies on its type III secretion system (T3SS) to inject effector proteins into host cells, making T3SS proteins attractive targets for antivirulence therapeutics. In this study, we employed an integrated computational pipeline to model and assemble the C. trachomatis T3SS constituent proteins. Template-based modeling using crystallographic structures of homologs from other Gram-negative bacteria revealed a highly conserved structural architecture despite low sequence identity (18-46%). Stereochemical validation confirmed high model quality, with most T3SS proteins exhibiting favorable protein-protein interactions (PPIs). Since the activity of the T3SS complex relies on extensive PPIs, we targeted these PPIs as a promising approach to attenuate bacterial virulence. CdsN, which functions as an ATPase of the T3SS, is a hexamer of which we targeted the dimerization interface. Structure-based virtual screening of compounds from the e-Drug3D and IMPPAT libraries against predicted hotspot residues and the identified druggable pocket at the CdsN dimeric interface, followed by ADMET screening, yielded three promising candidates: M Roflumilast (Drug ID: 1537), Elacestrant (Drug ID: 2081), and Tecovirimat (Drug ID: 1889). All three ligands formed thermodynamically stable complexes with the CdsN dimer, with Elacestrant demonstrating the most favourable binding free energy. This was also confirmed by 100 ns molecular dynamics simulation. This study provides new insights into the molecular architecture of C. trachomatis T3SS and identifies M Roflumilast, Elacestrant, and Tecovirimat as potential drug candidates against chlamydial infection. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=129 SRC="FIGDIR/small/723290v1_ufig1.gif" ALT="Figure 1"> View larger version (58K): org.highwire.dtl.DTLVardef@1821599org.highwire.dtl.DTLVardef@1581baaorg.highwire.dtl.DTLVardef@1805e98org.highwire.dtl.DTLVardef@c25e56_HPS_FORMAT_FIGEXP M_FIG C_FIG